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Abstract 



r-| ■ Hadronization corrections to the predictions of perturbative QCD are reviewed. 

J^. The existing models for the conversion of quarks and gluons into hadrons are 

Qh! summarized. The most successful models give a good description of the data 

(-H I on e"'"e~ event shapes and jet fragmentation functions, and suggest that the 

dominant hadronization effects have a 1/Q dependence on the hard process 
energy scale Q. In several cases the 1/Q terms can be understood in terms 
of a simple longitudinal phase-space model. They can also be inferred by re- 
5^ I lating non-perturbative renormalon effects to the infrared cutoff dependence of 

perturbative contributions. 
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1 Introduction 



Hadronic jet production at high energy coUiders has proved to be one of the most valuable 
testing-grounds for quantum chromodynamics (QCD). At high energies, or more precisely at 
large momentum transfers, the QCD coupling as becomes small and perturbation theory 
becomes more reliable. Perturbative predictions to next-to-leading order, and to higher 
order in a few cases, give a good account of the broad features of jet production processes 
in e"'"e~ hadron-hadron, and lepton-hadron collisions [|l],0. 

A serious barrier to further progress in jet physics, however, is our lack of understanding 
of the process of hadronization, in which the quarks and gluons of perturbative QCD are 
converted into the hadrons that are seen in the detectors. On general grounds, we expect 
that hadronization and other non-perturbative effects should give rise to power-suppressed 
corrections to quantities that are computable in perturbation theory. By power-suppressed, 
we mean corrections proportional to 1/Q^, where Q is the hard process scale (the centre- 
of-mass energy in e"'"e~ annihilation). At present there are no solid arguments to exclude 
contributions with p < 2 for observables like e"'"e~ event shapes, which are not fully inclusive 
with respect to final-state hadrons. Indeed, as we shall see, there are strong indications that 
the leading corrections are proportional to 1/Q in these cases. This is in contrast to deep 
inelastic structure functions and the total e"'"e~ hadronic cross section, where arguments 
based on the operator product expansion suggest that the dominant power corrections 
should decrease like l/Q"^ and l/Q^ respectively 0. 

The general picture of jet production and fragmentation which has developed and proved 
highly successful over the past decade has three fairly distinct stages [|]. First, on a 
scale of energy and time characterized by a large momentum transfer-squared Q^, a hard 
subprocess takes place involving a small number of primary partons (quarks and/or gluons). 
For example, in e"'"e~ annihilation the hard subprocess would be e"'"e~ -^ qq or, more 
rarely, e"'"e~ — >• qqg. Next, over a period characterized by scales t such that Q"^ > t > 
to, the primary partons develop into multi-parton cascades or showers by multiple gluon 
bremsstrahlung (Fig. |l|). These cascades, which tend to develop along the directions of 
the primary partons owing to the coUinear enhancements in QCD matrix elements, are 
the precursors of the jets that are observed experimentally. The showering cutoff scale, 
to, should be much greater than the intrinsic QCD scale A^, but is otherwise somewhat 
arbitrary, being set by the requirement that a perturbative description in terms of partons 
should remain appropriate down to that scale. 

After the parton shower has terminated, we are left with a set of partons with virtualities 
(virtual masses-squared) of the order of the cutoff scale to. From this point we enter the 
low momentum-transfer, long-distance regime in which non-perturbative effects cannot be 
neglected. The most important of these is hadronization, which converts the partons into 
the observed hadrons. As long as this stage of the process involves only small momentum 
transfers, presumably on a scale set by the QCD scale A ~ 250 MeV, it can be expected 
to lead to power corrections to quantities that are finite in perturbation theory, which are 
determined mainly by the two earlier stages of jet production and development. 

At present the only detailed descriptions of the hadronization process are provided by 
models, which are discussed briefly in the following section. In Sect. §, we consider the 




Figure 1: Parton cascade in e^e annihilation. 

predictions of these models for hadronization corrections in e"'"e~ event shapes, where a 
rather clear pattern of 1/Q corrections emerges. 

Next, in Sect. ^, we turn to a discussion of hadron energy spectra in jet fragmentation. 
Here hadronization enters in two distinct ways. First, the the shapes and normalizations 
of the spectra of the various hadronic species provide detailed tests of hadronization mod- 
els. Such quantities cannot be computed from perturbation theory, because they involve 
the observation of individual hadrons. However, once a fragmentation spectrum has been 
measured at a particular hard process energy scale Q, then its form at any other large 
scale can be predicted perturbatively using the factorization properties of QCD. The only 
ambiguity in the prediction, apart from the effect of higher-order terms that can be com- 
puted in principle, is again due to unknown power-suppressed corrections. We review the 
evidence on these corrections, which are poorly understood compared with the analogous 
(higher- twist) terms in deep inelastic scattering. 

One formal theoretical approach to power corrections that may be useful is the study 
of infrared renormalons, which are generated by the divergence of perturbation theory at 
high orders. Renormalons are associated with power-suppressed effects but it is not so clear 
what they have to do with the hadronization process. However, it appears likely that there 
is a renormalon contribution to event shapes, corresponding to a 1/Q correction. This is 
discussed in Sect. ^. Finally, in Sect. P, some conclusions are drawn. 



2 Hadronization models 

One general approach to hadronization, based on the observation that perturbation theory 
seems to work well down to rather low scales, is the hypothesis of local parton-hadron 
duality . Here one supposes only that the flow of momentum and quantum numbers at 
the hadron level tends to follow the flow established at the parton level. Thus, for example, 
the flavour of the quark initiating a jet should be found in a hadron near the jet axis. 
The extent to which the hadron flow deviates from the parton flow reflects the irreducible 
smearing of order A due to hadron formation. Perhaps the most striking example of local 



parton-hadron duality is the shape of the hadron spectrum in jet fragmentation at relatively 
low energies, which will be discussed in Sect. ^.4. 

The simplest more explicit hadronization model [^ is the longitudinal phase-space or 
'tube' model,|] in which a parton (or, more realistically, a colour-connected pair of partons) 
produces a jet of light hadrons which occupy a tube in (y,pj)-space, where y = ^ log[(-E + 
Pz)/{E — pz)] is rapidity and pt is transverse momentum, both measured with respect to 
the direction of the initial parton. If the hadron density in this space is p{pt), the energy 
and momentum of a tube of rapidity length Y are 



f 
E = dy d^ptp{pt)ptCosh.y = XsinhY 

Jo 

rY 

P = dy d'^ptp{pt)pt sinhy = AfcoshF — 1) ~ i? — A , (1) 

where A = / d'^ptp{pt)pt sets the hadronization scale. Notice that the jet momentum P 
receives a negative hadronization correction of relative order \/ E = 2X/Q for a two-jet 
configuration of total energy Q. Thus one generally expects hadronization effects to scale 
with energy like 1/Q. 

From Eqs. (|l]) we expect a mean-square hadronization contribution to jet masses of 

{m\^ = E^-P^^XQ. (2) 

Comparing the perturbative predictions for jet masses with experiment, one finds that a 
hadronization correction corresponding to 

A ~ 0.5 GeV (3) 

is required. Note that this implies a fairly large jet mass, about 7 GeV at Q ~ Mz, in 
addition to the perturbative contribution. 

We shall see that the above simple model successfully describes many of the gross 
features of hadronization. In order to make more detailed predictions, we need a specific 
model for the mechanism of hadronization. Over the years, three classes of models have 
been developed, which we outline briefly in the following subsections. 



2.1 Independent fragmentation model 

The simplest scheme for generating hadron distributions from those of partons is to suppose 
that each parton fragments independently. The original approach of Field and Feynman 
was designed to reproduce the limited transverse momenta and approximate scaling 
of energy fraction distributions observed in quark jets produced in e"'"e~ annihilation at 
moderate energies. The fragmenting quark is combined with an antiquark from a qq pair 
created out of the vacuum, to give a "first-generation" meson with energy fraction z. The 
leftover quark, with energy fraction {1 — z), is fragmented in the same way, and so on until 



TThis model is essentially the simplest version of the string model; we call it a tube to avoid confusion 
with the more sophisticated Lund string model discussed below. 



the leftover energy falls below some cutoff. Scaling follows from the energy independence of 
the distribution assumed for z, which is known as the fragmentation function. The limited 
transverse momenta come from the relative transverse momenta of the created qq pairs, 
which are given a Gaussian distribution. 

For gluon fragmentation, the gluon is first split into a quark-antiquark pair, either 
assigning all the gluon's momentum to one or the other (z = or 1) with equal probability 
P], so that the gluon behaves at a quark of random flavour, or using the g —>■ qq Altarelli- 
Parisi splitting function p. 

With about four parameters to describe the fragmentation function, the width of the 
transverse momentum distribution, the ratio of strange to nonstrange pair creation, and 
the ratio of vector to pseudoscalar meson production, the model proved quite successful in 
describing the broad features of two-jet and three-jet final states in e^e~ annihilation at 
moderate energies |8|-|T0|. 



A weakness of the independent fragmentation scheme, as formulated above, is that the 
fragmentation of a parton is supposed to depend on its energy rather than its virtuality. 
Indeed, the fragmenting parton is usually assumed to remain on mass shell, leading to 
violations of momentum conservation that have to be corrected by rescaling momenta after 
hadronization is completed. The residual colour and flavour of the leftover parton in each 
jet also have to be neutralized at this stage. There are further problems when jets become 
close together in angle. Instead of merging smoothly together into a single jet, as would 
happen if their fragmentation depended on their combined effective mass, even two precisely 
coUinear jets remain distinguishable from a single jet. 



2.2 String model 

The string model of hadronization |jll]-|14| is most easily described for e"'"e~ annihilation. 
Neglecting for the moment the possibility of gluon bremsstrahlung, the produced quark 
and antiquark move out in opposite directions, losing energy to the colour field, which 
is supposed to collapse into a stringlike configuration between them. The string has a 
uniform energy per unit length, corresponding to a linear quark confining potential, which 
is consistent with quarkonium spectroscopy. The string breaks up into hadron-sized pieces 
through spontaneous qq pair production in its intense colour field. 

In practice, the string fragmentation approach does not look very different from inde- 
pendent fragmentation for the simple quark-antiquark system. The string may be broken 
up starting at either the quark or the antiquark end, or both simultaneously (the breaking 
points have spacelike separations, so their temporal sequence is frame dependent), and it 
proceeds iteratively by qq pair creation, as in independent fragmentation. What one gains 
is a more consistent and covariant picture, together with some constraints on the frag- 



mentation function [13|, to ensure independence of whether one starts at the quark or the 



antiquark, and on the transverse momentum distribution |14], which is now related to the 



tunnelling mechanism by which qq pairs are created in the colour field of the string. 

The string model becomes more distinct from independent fragmentation when gluons 



are present [|T5[. These are supposed to produce kinks on the string, each initially carrying 
localized energy and momentum equal to that of its parent gluon. The fragmentation of 
the kinked string leads to an angular distribution of hadrons in e+e~ three-jet final states 
that is different from that predicted by independent fragmentation and in better agreement 
with experiment [|l^]. 

For multiparton states, such as those produced by parton showering at high Q^ (Fig. m), 
there is ambiguity about how strings should be connected between the various endpoints 
(quarks and antiquarks) and kinks (gluons). However, to leading order in N'"^ where 
A'^c = 3 is the number of colours, it is always possible to arrange the produced partons 
in a planar configuration, such that each has an equal and opposite colour to that of a 
neighbour (both neighbours, in the case of a gluon), like the quark and antiquark in the 
simplest e^e~ final state. The natural prescription is then to stretch the string between 
colour-connected neighbours, so as to make colour singlet strings of minimum invariant 
mass. The reformulation of parton showers in terms of sequential splitting of colour dipoles 



i r?! leads to the same rule for string connection. 



2.3 Cluster model 

An important property of the parton showering process is the preconfinement of colour 
||18|| . Preconfinement implies that the pairs of colour-connected neighbouring parton dis- 
cussed above have an asymptotic mass distribution that falls rapidly at high masses and is 
asymptotically Q^-independent and universal. This suggests a class of cluster hadroniza- 
tion models, in which colour-singlet clusters of partons form after the perturbative phase 
of jet development and then decay into the observed hadrons. 

The simplest way for colour-singlet clusters to form after the parton cascade is through 
non-perturbative splitting of gluons into qq pairs [113. Neighbouring colour-connected 



quarks and antiquarks can then combine into singlet clusters. The resulting cluster mass 
spectrum is again universal and steeply falling at high masses. Its precise form is determined 
by the QCD scale A, the parton shower cutoff to, and to a lesser extent the gluon-splitting 
mechanism. Typical cluster masses are normally two or three times \^. 

If a low value of the cutoff to is used, of the order of 1 GeV^ or less, most clusters 
have masses of up to a few GeV/c^ and it is reasonable to treat them as superpositions 
of meson resonances. In a popular model ||T9| , [20[1 , each such cluster is assumed to decay 



isotropically in its rest frame into a pair of hadrons, with branching ratios determined 
simply by density of states. The reduced phase space for cluster decay into heavy mesons 
and baryons is then sufficient to account fairly well for the multiplicities of the various 
kinds of hadrons observed in e"'"e~ final states. Furthermore, the hadronic energy and 
transverse momentum distributions agree well with experiment, without the introduction 
of any adjustable fragmentation functions. Also, the angular distribution in e^e~ three-jet 
events is successfully described, as in the string model [|16|, provided soft gluon coherence 



is taken into account in the parton shower via ordering of the opening angles of successive 
branchings pTl - p3| . 



An alternative approach to cluster formation and decay is to use a higher value of the 



cutoff to and an anisotropic, multihadron decay scheme for the resuhing heavy clusters [24]. 
Clearly, this approach lies somewhere between the low-mass cluster and string models. In 
practice, even with a low value of to one needs to invoke some such decay scheme for the 
small fraction of clusters that have masses of more that an few GeV/c^, for which isotropic 
two-body decay is an implausible hypothesis. 



2.4 QCD event generators 

There are three classes of programs for generating full events from parton showers, using 
each of the three above-mentioned hadronization models. 

Because of the difficulties discussed in connection with the independent fragmentation 
model, one would expect it to work best for the higher-momentum hadrons in final states 
consisting of a few well-separated jets. Probably on account of its simplicity, it was initially 
the model used most widely in conjunction with initial- and final-state parton showering 
for the simulation of hard hadron-hadron collisions, in the programs ISAJET [^ , CO JETS 



and FIELDAJET ^. 



The string hadronization model outlined above, with many further refinements, is the 



basis of the JETSET simulation program ||2^ , which also includes final-state parton show- 



ering with optional angular ordering. This program gives a very good detailed description 
of hadronic final states in e"'"e~ annihilation up to the highest energies studied so far [^ . 

The JETSET hadronization scheme is also used, in combination with initial- and final- 



state parton cascades, in the other very successful Lund simulation programs PYTHIA [30 



for hadron-hadron and LEPTO []5T| for lepton-hadron collisions. The alternative formula- 
tion of parton showers in terms of colour dipole splitting mentioned above is implemented 
in the program ARIADNE [P2[, which also uses JETSET for hadronization. 



The program HERWIG [^] uses a low-mass cluster hadronization model ^^ in con- 
junction with initial- and final-state parton cascades to simulate a wide variety of hard 
scattering processes. The showering algorithm includes angular ordering and azimuthal 
correlations due to coherence and gluon polarization. This approach gives a good account 
of diverse data with relatively few adjustable parameters [p9|. 



More detailed descriptions and comparisons of event generators for e"'"e physics may 
be found in Ref. [13^ 



3 Event shapes 

A popular way to study the jet-like characteristics of hadronic final states in e"'"e~ annihila- 
tion is to use event shape variables. The procedure is to define a quantity X which measures 
some particular aspect of the shape of the hadronic final states, for example whether the 
distribution of hadrons is pencil-like, planar, spherical etc. The differential cross section 
da/dX can be measured and compared with the theoretical prediction. For the latter to 



be calculable in perturbation theory, the variable should be infrared safe, i.e. insensitive 
to the emission of soft or coUinear particles. This is because the QCD matrix elements 
have singularities whenever a soft gluon or coUinear pair of massless partons is emitted. In 
particular, if p^ is any 3-momentum occurring in its definition, X must be invariant under 
the branching 

Pi -^ Pj + Pk (4) 

whenever p and pp. are parallel or one of them goes to zero. Quantities made out of linear 



sums of momenta meet this requirement. The most widely-used example is the thrust ||3^ 



^ = max , (5) 

2-ii \Pi\ 

where p^ are the final-state hadron (or parton) momenta and n is an arbitrary unit vector. 
If the p^ form an almost coUinear pair of jets, then after the specified maximization or 
minimization n will lie along the jets, defining the thrust axis of the event. As the jets 
become more pencil-like, the thrust approaches unity. 

A event shape variable that does not require finding an event axis is the C-parameter 

m 



In the case of a pencil-like two-jet event, the C-parameter is close to zero; the normalization 
is such that the maximum value of C is unity. 

As mentioned in Sect. |[4, the QCD Monte Carlo event generator programs HERWIG 
and JETSET are quite successful in describing the properties of e"'"e~ final states, including 
the distributions of event shape variables [^. From these successes, we can hope that 
the models built into those programs provide some guidance on the broad features of the 
hadronization process. What the programs suggest is that corrections to event shapes 
are still substantial at energy scales Q ~ Mz, typically around 10% of the leading-order 
QCD predictions. This is comparable with the next-to-leading 0(as) terms. Furthermore 
hadronization effects fall off rather slowly with increasing energy, apparently like 1/Q. 
Thus it is imperative to understand hadronization better in order to reap the benefit of 
any future 0{al) calculations of event shapes. 

We saw in Sect. ^ that 1/Q corrections are in fact generated by the simplest 'tube' 
model of hadronization. Consider for example the thrust distribution. The thrust of a 
two-jet event is precisely the jet momentum divided by the jet energy, and therefore in the 
tube model we expect a hadronization correction of —2X/Q. Thus for A ~ 500 MeV the 
hadronization correction to the thrust is expected to be 

The purely perturbative prediction for the mean thrust is ||38| 

(1 - T) = 0.335as + 1.02a^ + 0{al) , (8) 

which implies that (1 — T) ~ 0.055 at Q = Mz, assuming as{Mz) ~ 0.12. In fact, the value 
measured at LEP is (1 — T) = 0.068±0.003, consistent with an additional non-perturbative 
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Figure 2: Mean value of 1 - thrust in e"'"e annihilation. 



contribution of (1 GeV)/Q, which also agrees with the energy dependence of (1 — T) down 
to about Q = 15 GeV (Fig. 1). 

For the mean value of the C-parameter (j^), the tube model predicts a hadronization 
correction of SttX/Q, giving 

5 GeV 



{6C) 



had 



Q 



The perturbative prediction is 



{C) = 1.374as 



4.00a^ 



0{al] 



(9) 



(10) 



which implies (C) ~ 0.222 a.t Q = Mz, again assuming as{Mz) ~ 0.12. The value 
measured at LEP is (C) = 0.260 ± 0.006, which is consistent with the expected additional 
non-perturbative contribution, although in this case there are no published lower-energy 
measurements with which to check the Q-dependence. 

Turning to the differential distributions of event shape variables, hadronization seems 
mainly to cause a smearing by an amount proportional to 1/Q. The effect is therefore 
most pronounced in the neighbourhood of the peak of the distribution, as shown for the 
thrust distribution in Fig. |^. The shaded band shows the variation in the hadronization 
corrections deduced from JETSET and HERWIG, where the correction is defined as the 
ratio of the parton- and hadron-level Monte Carlo predictions. The predominant effect of 
hadronization is to smear out the parton-level peak at low values of 1 — T, enhancing the 
hadron-level distribution at intermediate values of the thrust. 
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Figure 3: Thrust distribution in Z° -^ hadrons, with detector and hadronization corrections. 

4 Jet fragmentation 



We have illustrated in Sect. |^ the use of perturbation theory to calculate 'infrared safe' 
quantities such as e"'"e~ event shapes. There are in addition predictions that can be made 
concerning some quantities that are 'infrared sensitive', i.e. that have infrared and collinear 
singularities in perturbation theory. Such quantities can still be handled provided the 
singularities can be collected into an overall factor which describes the sensitivity of the 
quantity to long-distance physics. The divergence of this factor corresponds to the fact 
that long-distance phenomena are not reliably predicted by perturbation theory. Therefore 
the divergent factor must be replaced by a finite factor determined either by experiment 
or according to some non-perturbative model. Once this is done, perturbation theory can 
be used to predict the scale dependence of the quantity. 

The best known example of a perturbative prediction concerning a factorizable infrared- 
sensitive quantity is the phenomenon of scaling violation in hadron structure functions. 
There one studies the parton distributions inside a hadron, probed by deep inelastic lep- 
ton scattering. Here we consider the related phenomenon for the fragmentation of a jet, 
produced for example in e"'"e~ annihilation, into hadrons of a given type h. 



4.1 Fragmentation functions 

The total fragmentation function for hadrons of type h in e"'"e~ annihilation at cm. energy 
^/s as is defined as 

F'^(x, s) = — ^(e+e- ^ hX) (11) 

(Ttot dx 

where x = 2Eh/\^ < 1 is the scaled hadron energy.| These functions are predicted by 
the QCD event generators discussed in Sect. |.4, and experimental data for the various 
identifiable hadron species provide strong constraints on the hadronization models used 
in the programs. Generally speaking, the most developed event generators, JETSET and 
HERWIG, give fairly good agreement with experiment, although there are problems with 
the yields of heavy strange particles p9|. 



The fragmentation function (|TT|) can be represented as a sum of contributions from the 
different primary partons i = u,d, . . . ,g: 

F\x, 5) = E / -C^i^■^ ^' c^s)D^{x/z, s) . (12) 

. J z 

In lowest order the coefficient function Cg for gluons is zero, while for quarks Cj = gi{s)S{l — 
z) where gi{s) is the appropriate electroweak coupling. 

We cannot compute the parton — > hadron fragmentation functions D^ from perturbation 
theory, since the production of hadrons is not a perturbative process. We might consider 
trying to compute functions Dl that would describe the fragmentation of partons of type i 
into partons of type j. However, they would be infinite, since the probability of emitting a 
collinear gluon or light quark is divergent. Nevertheless, it can be shown |^ that all such 
divergences are factorizable in the sense that one can write 

D{{x,t) = J2 f' dzK^iz,t,to) Diix/z,to) (13) 

k ■''' 

where the kernel function K^ is perturbatively calculable and K^{z,t,t) = ^^^^(l ~ z). 

In the real world, hadrons are formed and fragmentation functions are not divergent. As 
long as the scales t and to are large, we would not expect the form of Eq. (|T3|) to be affected 
by such long-distance phenomena. Therefore we replace j by h and apply the equation 
to hadronic fragmentation functions. After measuring these functions at some scale to, we 
can use the equation to predict their form at any other scale t, as long as both t and to are 
large enough for perturbation theory to be applicable in the calculation of the kernel K^. 



4.2 Scaling violation in jet fragmentation 

Consider the change in a fragmentation function Di{x,t) when the hard process scale is 
increased from t to t + dt. Such a change can only occur via the splitting of a parton of type 



■l-In practice, the approximation x = Xp — 2ph/\/s is often used. 
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i in this scale interval. Hence the fragmentation functions Di satisfy evolution equations 

(14) 



of the Altarelli-Parisi type, first introduced for deep inelastic scattering [^ 



t—Di{x,t) =J2J —^Pjiiz^^s)Djix/z,t) . 



The i —* j parton splitting function Pji has a perturbative expansion of the form 



,(0) 



"Sp(l), 



P,,{z,as) = P)':\z) + ^P)l\z) + 



27r 



(15) 



where the lowest-order splitting function P-j (z) is the same in fragmentation and deep 
inelastic scattering but the higher-order terms are different. The effect of splitting is qual- 
itatively the same in both cases: as the scale increases, one observes a scaling violation in 
which the x distribution is shifted towards lower values. 

The most common strategy for solving the evolution equations (0) is to take moments 
(Mellin transforms) with respect to x: 



D{N,t) 
the inverse Mellin transform being 



dxx'^~^ D{x,t) 



I 



D{x, t) = ^ I dN x-^ D{N, t) 
2-Ki Jc 



(16) 



(17) 



where the integration contour C in the complex A^ plane is parallel to the imaginary axis 
and to the right of all singularities of the integrand. After Mellin transformation, the 
convolution on the right-hand side of Eq. (|I^) becomes simply a product. 

The moments Pji of the splitting functions are called anomalous dimensions, usually 
denoted by 'yji{N, as). In lowest order they take the form 



7j?(iV) 



i'JiN) 



l^SiN) 



TSHiV) 



Cf 
Tr 
Cf 
2Ca 



1 1 

^ 2 ^ N{N + 1) 

{2 + N + N'^) 
N{N + 1){N + 2)^ 
'i2 + N + N'^) 



N ^ 
k=2 '^J 



N{m - 1) 

1 1 

^12 ^ A^(A^-l) 



;i8) 



Af 



t vi 

{N + l){N + 2) ^^k^ 



-:NfT, 



fJ-R , 



where the QCD colour factors are Cf = |, Ca = 3, Tr = i, and Nf is the number of quark 
flavours with masses less than the relevant scale \/i. 

We can consider fragmentation function combinations which are non-singlet in flavour 

Dg . In these combinations the mixing with the 



space, such as Dy = D^ 

flavour singlet gluons drops out and for a fixed value of as the solution is simply 



Da or Da 

H% HI 



Dv{N,t)=Dv{N,to){-J 



(19) 
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For a running coupling as{t), the scaling violation is no longer power-behaved in t. Inserting 
the lowest-order form for the running coupling, 

"^w = mm (20) 

where b = (IICa — 2Nf)/12TT, we find the solution 

Dy{N,t) = Dy{N,to) l^^j , dUN) = ^^ , (21) 

which varies like a power of In t. 

For the singlet quark fragmentation function 

Ds = j:iD,. + D,^), (22) 

i 

we have mixing with the fragmentation of the gluon and the evolution equation becomes a 
matrix relation of the form 

The anomalous dimension matrix in this equation has two real eigenvalues 7± given by 




1± = \ [igg + Iqq ± \J{lgg ' Iqqf + 8iV/7gg7gg] . (24) 

Expressing Ds and Dg as linear combinations of the corresponding eigenvectors i5+ and 
D_, we find that they evolve as superpositions of terms of the form (|21| ) with 7_|_ and 7_ 
in the place of 7gg. 

At small X, corresponding to A^ — *> 1, the g ^ g anomalous dimension becomes dominant 
and we find 7+ — > 7gg — > 00, 7_ -^ 7^^ — > 0. This region requires special treatment, as we 
discuss in the following two sections. 

There are several complications in the experimental study of scaling violation in jet 
fragmentation ||4^. First, the energy dependence of the electroweak couplings gi{s) that 
enter into Eq. (^) is especially strong in the energy region presently under study {^/s = 
20 — 90 GeV). In particular, the 6-quark contribution more than doubles in this range. The 
fragmentation of the b quark into charged hadrons, including the decay products of the b- 
flavoured hadron, is expected to be substantially softer than that of the other quarks, so its 
increased contribution can give rise to a 'fake' scaling violation that has nothing to do with 
QCD. A smaller, partially compensating effect is expected in charm fragmentation. These 
effects can be eliminated by extracting the b and c fragmentation functions at ^/s = Mi 
from tagged heavy quark events, and evolving them separately to other energies. 

Secondly, one requires the gluon fragmentation function Dg{x,s) in addition to those 
of the quarks. Although the gluon does not couple directly to the electroweak current, it 
contributes in higher order, and mixes with the quarks through evolution. Its fragmenta- 
tion can be studied in tagged heavy-quark three-jet {QQg) events, or via the longitudinal 
fragmentation function (see Sect. §.5). 

12 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0. 



Figure 4: Scaling violation in e"'"e fragmentation functions. 



A final complication is that power corrections to fragmentation functions, of the form 



f{x)/Q^, are not well understood. As in the case of event shapes, Monte Carlo studies PO 



suggest that hadronization can lead to 1/Q corrections. Therefore, possible contributions 
of this form should be included in the parametrization when fitting the scaling violation. 

Preliminary results of an analysis of scaling violation in the charged hadron spectrum by 
the ALEPH collaboration [^, based on a comparison of LEP data with those from lower- 
energy experiments, are shown in Fig. ^ Also included are their separate fragmentation 
functions for light (m, d, s) quarks, b quarks and gluons, showing that the latter two functions 
are significantly softer (quite similar to each other, within the errors) at a/s = Mz- An 
overall fit in the range 22 < ^/s < 91.2 GeV, 0.1 < a; < 0.8, incorporating full next- 
to- leading-order evolution and a simple parametrization of 1/Q power corrections, gives 
a good description of the data in the fitted region, as shown by the curves.^ The fitted 
power corrections are small and the value obtained for as(Mz) (0.127±0.011) is not highly 
sensitive to the form assumed for them. Indeed, a 1/Q^ dependence is not ruled out, and 
it may be that in this case the Monte Carlo models are misleading and 1/Q corrections 
are in fact absent. An operator-product approach to fragmentation, developed in Ref. [ ^ , 
does indeed suggest that 1/Q corrections are absent. 



§Thc fitted curves do not extrapolate well into the small- a: region, but a next-to-leading order treatment 
would not be expected to be reliable there: see Sect. 0.4. 
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4.3 Average multiplicity 

The average number of hadrons of type /i in a jet initiated by a parton i at scale t, Afl^{t), is 
just the integral of the fragmentation function, which is the N = 1 moment in the notation 
of Eq. (|1D: 

Mt{t)= f'dxD^{x,t) = b1{l,t). (25) 

JO 

If we try to compute the t dependence of this quantity, we immediately encounter the 
problem that the anomalous dimensions 7^^ and 7^3 given in Eq. (^) have poles at iV = 1. 
The reason is that for A^ < 1 the moment integrals are dominated by the region of small 
z, where Pgi{z) has a divergence associated with soft gluon emission. 

In fact, however, we can still solve the evolution equation for the average multiplicity 
provided we take into account the suppression of soft gluon emission due to QCD coher- 
ence p2|j23[| . The leading effect of coherence is obtained simply by changing the evolution 
variable from virtual mass-squared to the quantity t = E'^[l — cos^], where 6 is the opening 
angle, and imposing angular ordering, 9' < 9 [^,^. Thus in terms of the newly-defined 



evolution variable the ordering condition is t' < z t and Eq. (jT^) becomes 

t^A(x, t) = J2 [ v£^.-^(^' «s)I^,(x/z, zH) . (26) 

Notice that this differs from the conventional Altarelli-Parisi equation only in the z- 
dependent change of scale on the right-hand side. This change is not important for most 
values of x but we shall see that it is crucial at small x. 



For simplicity, consider first the solution of Eq. ( PB| ) taking as fixed and neglecting the 
sum over different branchings. Then taking moments as usual we have 

t^MN, t) = ^ f dz z''-'P{z)D{N, zH) . (27) 

at 2-71 Jx 

Now if we try a solution of the form 

D(Ar,t) oct^(^'"s) (28) 

we find that the anomalous dimension 7(A^, as) must satisfy the implicit equation 

^(iV, as) = — t dz ^A^-i+27(A',-s)p(^) . (29) 

27r JO 

When A^ — 1 is not small, we may neglect the 2'y{N, as) in the exponent of Eq. (^) 
and then we obtain the usual explicit formula for the anomalous dimension. For A^ ~ 1, 
the region we are interested in, the z —>■ behaviour of P{z) dominates and we may as a 
first approximation keep only those terms that are singular at z = 0. The most important 
such term appears in the gluon-gluon splitting function: Pgg{z) -^ 2Ca/ z as z — > 0. Then 
Eq. ( P^ implies that near A^ = 1 

,„(N, a,) = ^-^_^_A___ , (30) 
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and hence 



lgg{N, as) 



1 

4 



'(iV-l)2 + ^^^-(iV-l) 



TX 



27r 



27r 4^ ^ 32VCAas 



{N~lf + - 



(31) 



Thus for A^ — > 1 the gluon-gluon anomalous dimension behaves hke the square root of 
ag. How can this behaviour emerge from perturbation theory, which deals only in integer 
powers of as? The answer is that at any fixed A^ 7^ 1 we can expand Eq. ( |3TD in a different 
way for sufficiently small as: 



lgg{N,as) 



7r(iV-l) 



[CaoisX 



\ TT 



') (AT -1)3 



+ 



(32) 



This series displays the terms that are most singular as A^ — *> 1 in each order. These terms 
have been resummed in the expression (pi]), allowing the perturbation series to be analyti- 
cally continued outside its circle of convergence |as| < (7r/8C^)|A^ — Ip. By definition, the 
behaviour outside this circle (in particular, for the average multiplicity, at A^ = 1) cannot 
be represented by the series any more, even though it is fully implied by it. 

At sufficiently small x, the A^ ^ 1 singularity of the gluon-gluon anomalous dimension 
dominates in all fragmentation functions. Thus we obtain the behaviour (pSj) with 7 = 
'jgg for the total fragmentation function defined in Eq. ([121) . To predict this behaviour 
quantitatively we need to take account of the running of as, which can be done in the same 
way as for the other moments. Writing Eq. (^Sf ) in the form 



F{N,t) r^D{N,t) ocexp 






(33) 



we have to replace 'jgg{N, as) in the integrand by 'jgg^N, as(t'))- We then write 

-fgg{N,as{t))—= '' das , 



(34) 



where /3(as) = — ^a^ + ■ ■ ■, and find 



F{N,t) 



exp - 




(AT-l) 



1 27r , 



(35) 



Jag=ag(t) 



In e"'"e~ annihilation the scale t (the upper limit on E'^ [1 — cos 9] for any branching) is 
of the order of the centre-of-mass energy-squared s, and so the average multiplicity of any 
hadronic species has the asymptotic behaviour 



{n{s))=F{l,s) 




2Ca , / s 



(36) 
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Figure 5: Average multiplicity of charged hadrons in e"'"e annihilation. 



If the next-to-leading singularities of the anomalous dimensions, i.e. the terms with one 
less power of 1/(A^ — 1) in Eq. (|32D , are also resummed, the expression (|36|) is multiplied 
by a power of as{s) [^. The resulting prediction, shown by the solid curve in Fig. ^ is in 
very good agreement with experiment |^ . 



4.4 Small- a? fragmentation 



The behaviour of F{N, s) away from A^ = 1 determines the form of small-x fragmentation 
functions via the inverse Mellin transformation (p!?]). Keeping the first three terms in the 
Taylor expansion of the exponent, as displayed in Eq. (^5]), gives a simple Gaussian function 
of N which transforms into a Gaussian in the variable ^ = ln(l/a:;): 



xF{x, s) oc exp 



2a' 



i^-Cp 



where the peak position is 



^P 



1 



Abas{s) 



Ins 



and the width of the distribution of ^ is 



27r 



a 



m^CAalis) 



oc (Ins 



(37) 



(3^ 



(39) 



Again, one can compute next-to-leading corrections to these predictions ||^, and the re- 
sults agree very well with the form and energy dependence of fragmentation functions at 
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Figure 6: Distribution of ^ = ln(l/x) in e"'"e annihilation, 
small X as measured in e"'"e~ annihilation (Fig. ^) |^. This provides support for the no- 



tion of local parton-hadron duality mentioned in Sect. 0. We assumed in Eq. (|33D that 
the iV-dependence of the Mellin-transformed fragmentation function is dominated by the 
perturbative part involving the anomalous dimension, and that the non-perturbative factor 
is relatively smooth. Smoothness in A^ transforms into locality in x. Although the violation 
of scaling ensures that this is valid asymptotically, it seems to be a good approximation 
already at energies as low as 14 GeV. 

The energy-dependence (|38D of the peak in the ^-distribution is also a striking illustra- 
tion of soft gluon coherence, which is the origin of the suppression of hadron production at 
small X. Of course, a decrease at very small x is expected on purely kinematical grounds, 
but this would occur at particle energies proportional to their masses, i.e. at x oc m/y/s 
and hence ^ ~ ^ In s. Thus if the suppression were purely kinematic the peak position ^p 
would vary twice as rapidly with energy, which is ruled out by the data (Fig. ^. 



4.5 Longitudinal fragmentation 



Recently, the ALEPH [^] and OPAL |^ collaborations at LEP have presented preliminary 
results of the first analyses of the longitudinal fragmentation function in e"'"e~ annihilation. 
This is defined in terms of the joint distribution in the energy fraction x and the angle 6 



between the observed hadron and the incoming electron beam |40 



1 ,^2 Q Q Q 

{l + cos^e)FT{x) + -sin'^eFi^{x) + -coseFA{x) , (40) 



(Jtot dx d cos ^8 4 4 

17 



Figure 7: Position of the peak in the ^-distribution. 

where Ft, -Fl and Fa are respectively the transverse, longitudinal and asymmetric frag- 
mentation functions. [] Like the total fragmentation function, F = Ft + Fj^, each of these 
functions can be represented as a convolution of the parton fragmentation functions Di 
with appropriate coefficient functions C^ ' ' |^9[ as in Eq. (|12[). In fact the transverse and 



longitudinal coefficient functions are related in such a way that 



Fl{x) 



as 

2tx' 



C, 



dz 

z 



Ft(z)+A 



X 



- 1 1 D,iz) 



+ 0(a 



(41) 



Thus the gluon fragmentation function Dg can be extracted to leading order from measure- 
ments of Ft and Fl . 



Fig. ^ shows the preliminary OPAL data ||4^ on the transverse and longitudinal frag- 
mentation functions for charged particles. Note from Eq. (^TJ) that Fl is 0{as) relative 
to Ft; it also falls more steeply with x. Therefore even with LEP statistics the errors are 
large for x > 0.3. However, one still obtains useful information on the gluon fragmentation 
function over the full range of x, as shown in Fig. ^ Because the relation (^) is known only 
to leading order, there is an ambiguity about the energy scale at which Dg is measured. 
Comparisons with JETSET predictions at various scales are shown. 

Similar, although systematically somewhat lower, preliminary results on the longitu- 
dinal fragmentation function have been obtained by the ALEPH collaboration [41|, who 



include this information in their fit to scaling violation as a further constraint on the gluon 
fragmentation function shown in Fig. |^. 

Summed over all particle types, the total fragmentation function satisfies the energy 
sum rule 

- f dxxF{x) = 1 . (42) 

^All these functions also depend on the cm. energy y/s, which we take to be fixed here. 
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Figure 8: Transverse and longitudinal fragmentation functions in Z° -^ hadrons. 
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Figure 9: Gluon fragmentation function extracted from -Ft,l- 
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Figure 10: Longitudinal fraction of the e"'"e hadronic cross section. 



Similarly the integrals 



- / dxxFT^i^{x) = 



O"tot 



(43) 



tell us the transverse and longitudinal fractions of the total cross section. The perturbative 
prediction is 

^ = ^ + 0(ai). (44) 

O'tot TT 

that is, the whole of the 0{as) correction to (Xtot conies from the longitudinal part. Surpris- 
ingly, the 0(ag) correction has not yet been calculated. Once this has been done, Eq. (|^) 
will provide an interesting new way of measuring as. 

The preliminary OPAL data point for ai^/atot is shown, together with JETSET Monte 
Carlo predictions, in Fig. |1^. It should be noted that neither of the JETSET parton- 
level predictions fully includes the 0(ag) contribution |^. We see, however, that the 
data lie well above the leading-order prediction (dashed), which suggests that higher-order 
and/or non-perturbative corrections are significant. An estimate of the latter is provided 
by the difference between the JETSET curves for hadrons and partons. This difference 
shows a clear 1/Q behaviour, with a coefficient of about 1 GeV, illustrating once again the 
characteristic energy dependence of hadronization effects. 

In fact, in the simple 'tube' model of hadronization discussed in Sect. ^, the correction 
to the longitudinal cross is found to be 




had 



ttA 
2Q 



0.8 GeV 



(45) 



which agrees well with the JETSET prediction shown in Fig. [1^. The correction arises from 
mixing between the transverse and longitudinal angular dependences in Eq. (|40D due to 
hadronization. The transverse cross section receives an equal and opposite correction, and 
so there is of course no 1/Q term in the total cross section. 
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Figure 11: Mean charged multiplicity in e"'"e and ep current jet fragmentation. 

4.6 Jet fragmentation in deep inelastic scattering 

Jet physics in ep collisions promises to be a rich field of study for the new HERA collider 
at DESY. One can for example compare the properties of jets in e"'"e~ and ep collisions, 
which is a good test of the factorizability of jet fragmentation functions discussed above. 

The appropriate comparison is between one hemisphere of an e^e^ event, defined with 
respect to the thrust axis, and the "current jet hemisphere" of deep inelastic scattering in 
the Breit frame of reference P,p2|. The Breit frame is the one in which the momentum 
transfer from the lepton, q^, has only a ^-component: g^ = (0, 0, 0, Q). In the parton model, 
in this frame the struck quark enters from the right with 2;-momentum p^ = —\Q and sees 
the exchanged virtual boson as a "brick wall", from which it simply rebounds with p^ = 
+^Q. Thus the right-hand hemisphere of the final state should look just like one hemisphere 
of an e^e~ annihilation event at E^ra = Q- The left-hand ("beam jet") hemisphere, on the 
other hand, is different from e"'"e~ because it contains the proton remnant, moving with 
momentum p^ = — (1 — x)Q/2x in this frame. 

Fig. ^ shows the preliminary results of such a comparison by the ZEUS collaboration 



5T| , in this case comparing the charged multiplicity in the Breit frame current hemisphere 
with that in e"*"e~ annihilation, already shown in Fig. ^ We see that the two sets of data 
do seem to follow the same curve. 

A more elaborate comparison is shown in Fig. ^. Here the position of the peak in 
the distribution of ln(l/a;p), where Xp = 2p/Q in the Breit frame, is compared with that 
found in e+e" annihilation (Fig. 0). We see again that the preliminary HERA data follow 
and extrapolate the e"'"e~ curve, in remarkable agreement with the theoretical predictions 
derived in Sect. @.4. 
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Figure 12: Position of the peak in ln(l/xp) in e^e and ep current jet fragmentation. 

5 Power corrections from perturbation theory? 



We have seen that model studies suggest that for many infrared safe quantities, such as 
e"'"e~ event shapes, hadronization corrections are proportional to 1/Q where Q is the centre- 
of-mass energy. This is in contrast to the e"'"e~ total cross section, which for massless 
quarks has a leading power correction of order 1/Q^. The smallness of non-perturbative 
effects in the total cross section and related quantities has in fact led to a preference for 
such quantities as a means of determining as, even though event shapes have a stronger 
perturbative dependence on as- 

In the case of the total cross section, we also have some understanding of the leading 



power correction ||5^ . It is believed to arise from the vacuum expectation value of the gluon 
condensate, (asG^), which is the relevant gauge-invariant operator of lowest dimension, 
giving a correction proportional to {a^G'^) jQ'^. From this viewpoint, we do not even know 
why the corrections to event shapes should be of order l/Q: there are no local operators 
of dimension one to which they could be related. 



A possibly fruitful way of discussing power corrections is in terms of renormalons |53- 



55| . These are singularities of the Borel transform of the all-orders perturbative expression 
for a quantity, generated by factorial growth of the perturbation series at high orders. The 
idea of Borel transformation is that for any perturbation series 



^(«s) = EPn«S 



n+1 



(46) 



we can define a more convergent series 



n=0 



P(6) 






(47) 
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If this series converges we can define P{as) by 

P{as) = / dbe-^/^s P{b) (48) 







even when the original series is divergent. Then the series is said to be Borel summable. 
Probably the QCD perturbation series is not fully Borel summable 0, but we may still 
hope to resum important sets of contributions in this way. 

The growth of QCD matrix elements in the soft region is thought to give rise to a 
singularity (an infrared renormalon) in the 6-plane at the position h = 2 /bo (6o = [HC^a ~ 
2Nf]/12iT, called b in earlier sections), generating a power correction proportional to 



exp 



&o«s(Q2), 






The existence of a renormalon is supposed to indicate that the full QCD prediction would 
exhibit the same power correction. In this language, the appearance of 1/Q corrections 
to event shapes would be associated with a new infrared renormalon at l/26o in the Borel 
transforms of these quantities. 



In a recent paper |5^, an idea about renormalons due to Bigi, Shifman, Uraltsev and 



Vainshtein ^E\ was applied to e+e" event shapes. Those authors suggested that there is a 
simple correspondence between infrared renormalon positions and the power corrections to 
fixed-order perturbative predictions evaluated with an infrared cutoff. In e"'"e~— >■ hadrons, 
to first order in ag all diagrams are QED-like and a suitable cutoff can be imposed by 
introducing a small mass /i (in principle much greater than A) into the gluon propagator.^] 
According to Ref. |5^, the first-order perturbative total cross section with such a cutoff, 
normalized to the Born value, takes the form 

Rp = l + as/7r-DasijyQ^ + ... (50) 

where D is a constant and the ellipsis represents non-leading power corrections. The 
dependence on fi must cancel between this expression and the soft contribution, which 
builds the renormalon. Thus the leading renormalon has to occur at 2/6o, as stated above, 
and we expect a non-perturbative contribution of the form 

i?„p=[CA^ + Z}as(/x)/i1/Q' (51) 

where C is a constant. The /i-dependence appears as an arbitrariness in the part of the 
correction that we attribute to the renormalon and the part that is generated in fixed order. 

In the case of event shapes, however, one finds that the introduction of a gluon mass 
leads to corrections of order ag/x/Q, which cancel in the total cross section. That is, for a 
generic (infrared safe) event shape S one has in first order 

Sp = Asas- Ds as ^-/Q + • • • (52) 



instead of a result of the form (|50D . The coefficients Ds are easily computed; to this order 
they arise entirely from the reduction of phase space for real gluon emission. The values 
obtained in Ref. [^| for some representative quantities are listed, together with the leading 
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s 


As 


Ds 


Cs^ 


(l-T) 


0.335 


i = l-7 


~ 1.0 GeV 


(C) 


1.375 


8 


~ 5.0 GeV 


(yi 


I/tt 


4/3 


~ 0.8 GeV 



Table 1: Coefficients of terms in Eqs. ([5^) and (p3| 



coefficients As-, in Table [1|. Here T is the thrust, C is the C-paranieter, and 0"^ is the 
longitudinal cross section, as defined in earher sections. 

From Eq. ( |5^ ) we expect that a new infrared renormalon at l/26o is present in event 
shapes, leading to a non-perturbative contribution 



Snp = [C5 A + Ds as(/i) yu]/Q , 



(53) 



whose dependence on the arbitrary cutoff /i cancels against that of the perturbative part, 
leaving a cutoff-independent power correction CsA/Q. The observed values of these cor- 
rections, given by Eqs. (|^), @ and (^5]), are also shown in Table 1. 

It is remarkable that the observed 1/Q corrections are, within the uncertainties, pro- 
portional to the perturbative coefficients Ds, suggesting that Eq. (|53D takes the general 
form 



5', 



np 






l + rfas(/i)-r 



(54) 



where d is a. constant. In fact, for the quantities shown, the cutoff dependence is propor- 
tional to the hadronization correction computed using the simple 'tube' model introduced 
in Sect. ||. Thus one could argue that the success of that model is due to some underlying 
connection with the renormalon contribution. Similarly the success of the Monte Carlo 
hadronization models could be a reflection of the fact that introducing an infrared cutoff 
on the parton shower (called to in Sect. H) is sufficient to reproduce the systematics of the 
leading power corrections to event shapes. 

It would be interesting to apply the above approach to a wide variety of quantities, and 
to try to construct event shapes from which 1/Q corrections are absent. From Table 1 we 
see that the combination T -|- 2C/37r might be of this type. It would be desirable to extend 
the treatment to higher orders in perturbation theory, but it is difficult to see how this can 
be done in a gauge-invariant wayf*] 



II Recall that in QED a photon mass can be introduced in this way without violating the Ward identities 
associated with current conservation; see Ref. [£^, p. 136. 

**The cutoff procedure used in the Monte Carlo models corresponds implicitly to the choice of a particular 
axial gauge. 
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6 Conclusions 



One firm conclusion that can be drawn about hadronization is that we are still far from 
understanding it. More experimental information is needed: it would be particularly valu- 
able to have data on power corrections for a wide range of e"'"e~ observables, along the 
lines of Fig. |^ for thrust, with proper analysis of errors on the fitted powers of Q and 
their coefficients. Analyses of hadronization corrections to event shape distributions in 
terms of smearing rather than correction factors would be useful: how do {^T)^^ and 

((ST)'^ — {6T) ) depend on T and Q? Some information on Q-dependence will be ob- 

\ / had 

tained from LEP2, but it would be sensible and profitable to go back and re-analyse as 
much lower-energy data as possible in new ways. 

The motivation for further experimental and theoretical study of hadronization is twofold. 
First, after a long period of inactivity, apart from some model building, there has recently 
been a revival of theoretical interest in power corrections to perturbative predictions in gen- 
eral. The interest stems mainly from new ideas about the behaviour of perturbation theory 
at high orders, and about the relationship between perturbative and non-perturbative ef- 
fects. Thus there are new theoretical conjectures to be tested experimentally. 

The second motivation is more pragmatic: 0{a^) calculations of more and more ob- 
servables are steadily becoming available. To use the extra power of these predictions to 
measure as we need better control of power corrections. In particular, it would be useful if 
1/Q corrections could be calculated, or if combinations of observables could be found from 
which they are absent. As an example, it was conjectured in Sect, ^that this might be the 
case for the quantity T + 2C/37r. 
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